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Abstract 

We compare experimental spectra of radiative and semileptonic B decays with the predictions of a model 
based on soft-gluon resummation to next-to-next-to-leading order and on a ghost-less time-like coupling. We 
find a good agreement with photon spectra in the radiative decay and with hadron mass distributions in the 
semileptonic one: the extracted values for as (vnz) are in agreement with the current PDG average within at 
most two standard deviations. The agreement is instead less good for the electron spectra measured by BaBar 
and Belle in semileptonic decays for small electron energies (< 2.2 GeV): our spectrum is harder. We also show 
that, in general, the inclusion of next-to-next-to-leading order effects is crucial for bringing the model closer to 
the data and that the non-power expansion introduced in the framework of analytic coupling studies does not 
accurately describe soft-gluon effects. 
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1 Introduction 



The aim of this work is to analyze measured B decay spectra with a model based on (i) soft gluon resummation 
to next-to-next-to- leading order and (ii) an effective QCD coupling having no Landau pole £Q. This coupling is 
constructed by means of an extrapolation at low energy of the high-energy behavior of the standard coupling. 
More technically, an analyticity principle is used. 

B decay spectra are substantially affected by long-distance effects, the most important ones being the soft 
interactions occurring in the fragmentation of the B meson into the b quark. The B meson — a colorless 
composite particle — emits the spectator quark and radiates soft gluons, 

B — ► b + sp + g\ + g 2 H + g n , (1) 

to convert into the colored b quark which later decays because of weak interactions, 

b -> s + 7 (2) 

or 

b -» u + I + v. (3) 

Perturbation theory can describe the fragmentation of a & quark into a b quark with a fraction of the original 
energy-momentum as an effect of multiple gluon radiation, but it clearly cannot describe that part of the 
fragmentation involving the spectator quark. Physical intuition suggests that initial bound-state effects are 
substantial for 

m\ « m B A QC D « 2GeV 2 , (4) 

which is experimentally interesting: that is the well known Fermi motion of the b quark in the B meson {mx 
is the final invariant hadron mass). This non-perturbative effect — which classically can be pictured as a small 
vibration of the b quark in the B meson because of the interactions with the spectator — has been formalized in 
an effective field theory by means of the well-known shape function or structure function of the heavy flavors [2]. 
Many models have been constructed to describe Fermi-motion as a genuinely non-perturbative effect involving 
the hadron structure perturbative corrections are included, if desired, later on and play in any case a 
minor role. In this work we adopt a different philosophy: in essence, we assume that the fragmentation of 
the lowest-lying beauty meson into the beauty quark and the spectator quark can be described as a radiation 
process off the b with a proper coupling. Even though dynamics of light degrees of freedom in the B meson 
is complicated, we assume that the related effects on semi-inclusive spectra are simple. More precisely, we 
assume that bound-state effects can be incorporated into an effective QCD coupling, which is inserted in the 
standard soft-gluon resummation formulas. We extrapolate therefore the perturbative QCD formulas to a non- 
perturbative region by assuming that the relevant non-perturbative effects can be relegated into an effective 
coupling. Since the perturbative formulas involve truncated expansions in the QCD coupling, it is clear that our 
approach is meaningful as long as the effective coupling remains appreciably smaller than one in all the relevant 
integration range. From Fig. 0we see that our effective coupling is « 0.5 for a typical soft scale k± ss 0.5 GcV 
(corresponding to x 7 = 2E 7 /ttib ~ 0.9 in radiative B decays), i.e. it is reasonably smaller than one. 

Since the whole fragmentation process is described in a perturbative framework, we do not distinguish 
between the mass of the B meson and the pole mass of the b quark, i.e. we consistently set rrib = tub- We also 
assume that this effective coupling is universal, i.e. that it can be used to describe different processes, and that 
it can be constructed on the basis of analyticity arguments. These are additional assumptions with respect to 
the basic one, which could eventually be relaxed. 

Let us remark that the resummed perturbative expansion for semi-inclusive quantities is incomplete even at 
the formal level. For inclusive quantities characterized by a hard scale Q, the cross section can be written in a 
consistent way as an expansion in the coupling at the scale Q, 

oo 

<WQ) = 5>„ag(Q 2 ), (5) 

71=0 
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where the c„'s are numerical coefficients of order one: no prescription is needed. Semi-inclusive processes are 
instead multi-scale processes, characterized by fluctuations with transverse momenta up to Q; the physical origin 
is very clear: a jet with a relatively large invariant mass mx (Aqcd *C rnx <C Q) can contain very soft partons, 
with transverse momenta of the order of the hadronic scale. Unlike case (JSJ, one has to face perturbative 
contributions of the form 

f Q2 rik 2 

/ (6) 

where an ill-defined integration over the Landau pole is made, even for large Q 3> Aqcd- A prescription for 
the low-energy behaviour of the coupling is therefore needed in any case. Even if quark confinement did not 
exist and partons instead of hadrons were the asymptotic states, a prescription would anyway be necessary to 
compute resummed cross-sections. 

It is clear that our approach has intrinsic and obvious limitations. The mass of the proton, for example, 
cannot clearly be computed by means of perturbative formulas with an effective coupling inserted in them: a 
genuinely non-perturbative technique is mandatory in this case, such as lattice QCD. Our point is that, with 
an effective coupling, we want to describe Fermi motion only, i.e. a specific non-perturbative effect, not all 
non-perturbative effects. We do not aim for example at describing the K* peak which appears in the radiative 
hadron mass distribution (see Fig.[3J), or, equivalently, the ir and p peaks which appear in the semileptonic one 
(see Fig.[7]and Fig. 0. These peaks, occurring for 

™ 2 x « Si 2 QCD, (7) 

are related to final-state hadronization, i.e. to the recombination of partons into hadrons. This effect has a 
different nature with respect to Fermi motion and occurs at a different scale (cfr. eq. (@J with eq. (JJJ)). With 
our model, we just want to describe a broad peak in the hadron mass distribution occurring in region ifljl. A 
possible difference between the photon spectra of, let's say, B — » X s ^f and A& — > Xg'y decays, could not be 
described or naturally incorporated in our model, which is a kind of "spectator model for spectra" . 

The validity of our approach cannot be judged a priori, but only a posteriori, by comparing its predictions 
with experimental data. One may ask which is the advantage of our approach compared to the standard 
one of postulating directly shapes for the non-perturbative components of the spectra and convolving them 
with the perturbative ones in the minimal prescription [I]. The answer is that we want to take advantage of 
the universality properties of QCD radiation, which are reflected in rcsummation formulas. In the standard 
approach, one has to postulate ad-hoc and un-related shapes for the non-perturbative components entering 
different observables, such as heavy flavor decay spectra, heavy flavor fragmentation, e + e _ shape variables, 
etc. If universal aspects of QCD dynamics — as measured in different processes — do exist, such aspects are 
not easily uncovered with the standard approach. On the other hand, with our method, such an investigation 
looks rather natural: to describe different processes, we use different perturbative formulas — quite often the 
same formulas but with different coefficients — with the same effective coupling by assumption and we look 
at the data jSJ- Our philosophy involves a "one step" approach: we deal simultaneously with perturbative and 
non-perturbative effects. The standard method is instead a "two-step" approach: one resums the perturbative 
long-distance effects in a minimal way — picking up just the infrared logarithms — and then introduces a 
physically motivated non-perturbative model. 

Another advantage of our approach is that it allows for a simple extraction of the value of the standard 
QCD coupling at a reference scale, f.i. as (mz), by comparing its predictions with measured B decay spectra. 
That is because the model uses ordinary perturbative formulas with a prescription for the coupling in the 
low-energy tail, and therefore there is not any double-counting problem in merging together short-distance and 
long-distance effects. A peculiarity of our model is that it has no free parameters, apart of course the true QCD 
ones, i.e. the hadronic scale Aqcd and the quark masses m q 's. It is therefore "rigid", in the sense that there 
is not a natural way to tune it to fit the data. 

The plan of the paper is as follows. 

In sec. El we summarize the main features of the ghost-less QCD coupling, which is basically an extrapolation 
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of the ordinary QCD coupling down to small momentum scales according to an analyticity principle which 
removes the Landau pole. 

In sec.|3]we construct the effective coupling controlling the evolution of gluon cascades, which are intrinsically 
time-like processes. The absorptive effects related to the decay of the time-like gluons are included in this 
effective coupling to all orders in perturbation theory. 

Sec.0]is the main one and describes the model based on soft gluon resummation in NNLO and on the effective 
coupling constructed in the previous section. A discussion of the relevance of the next-to-next-to-leading-order 
effects in our model is also presented. We also comment on the non-power expansion introduced in analytic 
coupling studies. 

In sec. we apply the model to describe B — > X s j decays. We compare its predictions with the invariant 
hadron mass distribution measured by BaBar and with the inclusive photon spectrum measured by Cleo, BaBar 
and Belle. Since these spectra are independent from each other, we obtain for each of them a value of as{mz) 
which optimizes the agreement with the data. 

In sec. we apply the model to the charmless semileptonic decays B — » X u lv. We compare our predictions 
with the invariant hadron mass distribution measured by BaBar and Belle and with the charged lepton energy 
spectrum measured by Cleo, BaBar and Belle. We extract values of as(mz) as discussed above. 

Finally, in sec. we draw our conclusions concerning the agreement of the model with the data. We also 
consider natural developments and improvements. 

There is also an appendix collecting formulas for the radiative decay and an appendix with tables of values 
of the QCD form factor in our model for a set of values of as(mz)- 



2 Ghost-less Coupling 

Let us begin considering QCD regularized with an ultra-violet cut-off Ao and with a bare coupling ao. The 
correlation function 4 representing the quark-gluon interaction has a perturbative expansion of the form: 

A 2 A 2 
T W9 (Pi = Pi = Pa = I 2 ) = a o + flo «o lQ g _ 2 —■ + a l c + Pi a l lQ - 2 ' " 



I -/3 a log A 2 /(-q*-ie)^ [> 

where for simplicity's sake we have considered the symmetric point p\ = v\ — v\- I n the- last member we 
have resummed the well-known geometrical series of the leading logarithms. /3q = (11 — 2/3 nf) /(4-7r) is the 
first-coefficient of the /3-function, nf is the number of active flavors and c is a real constant whose explicit 
expression is not relevant here. This Green function can be used to define the renormalized QCD coupling [HJ: 

r * (PI = P ^ p 3 = ^) * i-ftao log "III (-, 2 - ie) S a{ ~ q2) - (9) 

To have a real coupling, one generally assumes a space-like configuration of the momenta, 

q 2 < (10) 

and to avoid explicit minus signs in the renormalization conditions, one defines, like in Deep-Inelastic-Scattering 
(DIS): 5 

Q 2 = -Q 2 - (12) 

4 To be accurate, we consider the qqg correlation function amputated of all legs and written in terms of the renormalized fields. 
5 Note that Q 2 > in the space-like region while Q 2 < in the time-like one and the e-prescription for Q 2 is opposite to that 
for q 2 : 

Q 2 = Q 2 - ie. (11) 
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We then obtain the usual expression for the renormalized QCD coupling in leading order (LO): 

aio[Q2) = l-A)aologAg/Q2 = A)logQ 2 /A^ CD ' (13) 
where on the last member we have introduced the QCD scale 

1 



A QCD = A ex P 



(14) 



For notational simplicity, let us write A in place of Aqcd from now on. The function on the r.h.s. of eq. i|13fl 
has: 

1. a cut for Q 2 < 6 , related to the decay of a time- like gluon into secondary partons, 

g* -> 99, qq, ■■■ ■ (15) 

This singularity has therefore a clear physical meaning; 

2. a simple pole for Q 2 — A 2 , which does not have any physical meaning 7 . This singularity is often called 
"Landau ghost" because of its original appearance in QED in the interacting electron propagator |S] . It 
implies a formal divergence of the coupling and a breakdown of the perturbative scheme. 

It has been suggested to replace the usual expression for the coupling in eq. (| 1 31) with a "ghost-less" or "analytic" 
coupling a having the following properties pQ: 

1. it has the same discontinuity along the cut as the standard coupling: 

Disc a = Disc a (time — like region); (16) 

2. it is analytic elsewhere in the complex plane. 
Let us now consider the function 



/(-) - ^ (17) 



s + Q 2 

where Q 2 is a complex number not lying on the negative axis including the origin. By assumption, /(s) is 
analytic in the complex s-plane cut along the positive axis s > 0, except for a simple pole in 



s = -Q z . (18) 

We apply the residue theorem to f(s) integrated along a closed contour F avoiding the "physical" cut for s > 0, 
containing a circle of infinitesimal radius around the origin c e (e — ► 0), a circle at infinity c r (r — > oo), a line 
above the cut (s — > s + ie) and a line below the cut (s — ► s — is): see Fig. 2] Being the pole l(T8jl the only 
singularity inside the contour, we obtain the following expression for the analytic coupling: 7 



2m Jrs + Q 2 2iti J c s + Q 2 2m J c s + Q 2 



1 r ^- s - i£ y ds + ± f a{ - s+ *td S . d9) 



2m J s + ie + Q 2 2m s — ie + Q 



6 As usual, the logarithm function is cut along the negative semi-axis, so that: log(— 1 ± ie) = ±itt. 

7 Note that, had we taken Q 2 < 0, the pole 1181 would have been located on the cut and the integral of f(s) over F would have 
been zero. 
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Figure 1: integration contour T used to construct the ghost-less coupling. 



We assume that the contributions of c e and of c r vanish. Since s + Q 2 ^ for s > 0, 



lim 

S^0 + 



a(— s — ie) 



a(—s + ie) 



%e 



is 



1 



lim [a(— s — ie) — a{— s + ie)] 



_^Disc s a(- S ), 



S + Q 2 A"5+ 

where the discontinuity of a function F(s) is defined in general as: 

Disc s F(s) = lim \F(s + ie) - F(s-ie)]. 

e->0+ 

Taking into account that for s > (see eq. I|16|l) 

Disc s a(— s) = Disc s a(—s), 
we obtain the following integral representation for the ghost-less coupling in terms of the standard one: 

1 r ds 



2ni 



Q- 



Disc s a(—s). 



(20) 
(21) 

(22) 

(23) 



Eq. (|23|l is just a dispersion relation which, for clarity's sake, has been fully derived from first principles. By 
inserting on the last member the expression for the standard coupling at lowest order as given by eq. I|13[l . we 
obtain: 



aio(Q 2 ) = lim 



1 



ds 



o+ 2nil3o J s + Q 2 



1 



log(-s/A 2 - ie) log (-s/A 2 + ie) 



(24) 



The integral above is elementary. It can also be computed with the residue theorem by considering the contour T 
above. The circle of infinitesimal radius around the origin and the circle at infinity give vanishing contributions 
to the integral. There are two simple poles in s = — Q 2 and in s = — A 2 , so that: 



aUQ 2 ) = -r 



l 



A 2 



0o LlogQ 2 /A 2 Q 2 -A 2 



(25) 



Let us make a few remarks: 



1. comparing the r.h.s. of eqs. (|13|l and (|25|l . we see that the "analyticization" procedure had the effect of 
subtracting the infrared pole in Q 2 = A 2 by means of a power-suppressed term, in a minimal way; 
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2. the analytic coupling has a constant limit at zero momentum transfer: 



lim & lo (Q 2 ) = -j- « 0(1); (26) 



3. the term added to the standard coupling, 

1 A 2 



f3 Q 2 - A 2 ' 



(27) 



does not modify the high-energy behavior because it decays as an inverse power of the hard scale, i.e. 
infinitely faster than any inverse power of the logarithm of Q 2 . In more formal terms, the added term 
l|27() is exponentially small in the coupling, and therefore is always missed in an asymptotic expansion for 
Q 2 -> oo: 



A 1 -l/[/3 a !o (Q 2 )] 



Q2 _ ^2 e l/[/3oaio(Q 2 )] _ 1 

4. since the power correction has no discontinuity in the time-like region Q 2 < 0, 



(28) 



Disc Q2 A _ 2 A 2 = for Q 2 < °> ( 29 ) 

it trivially follows that the analytic coupling has the same discontinuity as the standard one, as originally 
requested. 



Let us now discuss the extension to next-to-leading order (NLO). The NLO correction to the standard coupling 

^>--W*£)i3W (30) 

where p\ is the second-order coefficient of the /3-function in the normalization assumed in [5], involves: 



1. the factor 1/ log 2 Q 2 /A 2 , having a cut for Q 2 < 0, related to the decay of the time-like gluon into on-shell 
partons (see eq. |TBJ), and a double pole for Q 2 = A 2 ; 

2. the factor log (logQ 2 /A 2 ), having a cut for Q 2 < related to the "internal" logarithm and another cut 
for < Q 2 < A 2 related to the "external" logarithm. 



The singularities for Q 2 — A 2 and for < Q 2 < A 2 are unphysical because they refer to the space-like region, 
where the virtual gluon cannot decay into physical parton states. "Analyticization" can be made as in lowest 
order: one requires that the analytic correction term has the same discontinuity for Q 2 < as the standard 
one but it is regular elsewhere in the complex plane: 

MO 2 )^-^ —Dtac.faO). (31) 

The integral above — unlike the lowest-order case — is not elementary but it can easily be made numerically. 
The following remarks are in order: 



• the value of the analytic coupling at zero momentum transfer is not modified in higher order because: 

lim 6a(Q 2 ) = 0; (32) 

• it can be shown that Sa(Q 2 ) has the same logarithmic terms as 5a(Q 2 ) ^H], i-e. that the difference resides 
in power-suppressed terms, as we have explicitly found for the leading order. 
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The NLO coupling is defined as: a n i — ai a + 5a. Within our accuracy, the next-to-next-to-leading order 
(NNLO) corrections to the coupling are also needed: 



log 2 I log I - log I log 



Q \ . L Q'\ , Pah-Pi 



A 2 J b V A 2 J pf 



1 



log 3 Q 2 /A 2 ' 



(33) 



where 02 is the third-order coefficient of the /3-function. One finds similar singularities as in the NLO case, 
which are removed again according to principle of "minimal analyticity" already used: 

5a' (Q 2 ) = — / -^Di S c s <W(- S ). (34) 

Jo s T W 

The NNLO analytic coupling reads (see Fig. |5J : 

ol = olio + Sa + 5a'. (35) 

Let us remark that an expansion in powers of the analytic coupling a is an asymptotic expansion — as in the 
standard case — because the logarithmic structure of a is the same as that of a. As Fig. [3 clearly shows, the 
standard coupling and the ghost-less one are barely distinguishable at large scales PQ. 



3 Effective Coupling for Gluon Cascade 

As well known from perturbation theory, the emission of a gluon in a process is accompanied by an additional 
factor a in the cross section, where a is the tree-level QCD coupling. In higher orders, one has to consider: 

1. multiple emissions off the primary color charges — the heavy and the light quark in B decays; 

2. secondary emissions off the radiated gluons. 

Primary multiple emissions produce the exponentiation of the one-gluon distribution while secondary emissions 
produce the decay of the radiated gluons into secondary partons — see eq. I)15[l. In the case of form factors, 
which are inclusive with respect to gluon decays, these higher-order terms have the main effect of replacing the 
tree-level coupling with an effective coupling evaluated at the transverse momentum of the primary emitted 
gluon [TT]: 

a -> a(k 2 ± ), (36) 
where 2 

a{kl) = I X dsDisc,^^. (37) 



2tt Jo 

The coupling a(k±) is characteristic of the QCD cascade and it is given by the integral of the discontinuity of 
the (interacting) gluon propagator over virtualities s cut-off by the primary gluon transverse momentum. Let 
us remark that the cascade (or effective or time-like) coupling always refers to time-like kinematics. 

The prescription at the root of our model is simply to replace the standard coupling on the r.h.s. of eq. I|37|l 
with the ghost-less coupling constructed in the previous section: 

a(kl) = 1- f k± dsDisc,^^. (38) 



2vr Jo 

If we neglect the — iir terms in the integral over the discontinuity — i.e. the absorptive effects — the cascade 
coupling exactly reduces to the ghost-less one: 

a(kl) a{k\). (39) 
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To render our model as accurate as possible, we include such absorptive effects and perform the integral on the 
r.h.s. of eq. exactly. By inserting the analytic coupling at LO in the integrand on the r.h.s. of eq. P5|l. we 
obtain for the effective coupling: 



1 



2m(3 

At NLO, one has to add the contribution: 



log log 



kl 

A 2 



5a(kl) = 



Ih 



2ni(3§ 



lOgfc^/A 2 + ITT 



logfc^/A 2 



log ( log k\ /A 2 + i-xj + 1 log ( log k\ /A 2 -mj+l 



(40) 



(41) 



The NNLO corrections read: 



6a'(k 2 x ) 



Pi 



log 2 (log k 2 ± /A 2 + in) log 2 (log kl I A 2 - in) 



(logfc 2 /A 2 + inY 



(log*; 2 /A 2 -mY 



/3 2 -/?o/3 2 



(log/c 2 /A 2 + z7r)' {\ogk 2 J A 2 - ittY 



The time-like coupling in NNLO is simply the sum of the above terms: 

a(k ± ) = ai (kl) + Sa(kl) + Sa'(kl). 

Let us make a few remarks: 



(42) 



(43) 



1. the cascade coupling is very close to the ghost-less one for very small scales, let's say less than 1 GeV 
(see Fig. |2J . That is partly a consequence of the fact that both couplings have the same limit at zero 
momentum, l//?ch and is partly accidental [J. The cascade coupling is instead smaller than the standard 
coupling or the ghost-less one in the perturbative region, at large scales, because it has an additional 
negative third-order contribution w — 1/ log 3 Q 2 — see next point; 



2. the time-like coupling a has an expansion in powers of the standard MS coupling a of the form: 



a 



5 j3 

6 A) 



- (7r/? ) 2 a 4 + 0{a 5 ). 



(44) 



The relation above can be considered as an ordinary change of scheme for the coupling starting at third 
order; 



3. the (3 function for the time-like coupling, defined by the relation 

da 



d\og kl 

has a negative third-order coefficient 8 



[3(a) 



f3 Q a 2 - 0i a 3 - /3 2 a 



^3 



^.4 



h = (h - I (7r/3o) 2 /3o < 0, 
in agreement with the fact that the coupling saturates at small scales; 



(45) 



(46) 



4. An expansion in powers of a(kl) is not an asymptotic expansion for kl — > oo because a even at LO 
contains infinitely many inverse powers of log kl . 

8 The first two coefficients are, as well known, invariant under a change of scheme: /3q = /3q, Pi = Pi- 
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The decoupling relations for the time-like coupling differ from the ones for the standard MS coupling and read: 



11 

72^ 



17 



54 54 



/ 1 Z.3 



(47) 



where a n , is the time-like coupling with n/ active flavors and a nf —i with n/ — 1. The above relation has to be 
imposed at a scale /j, such that m(/z) = fj,, where m(/i) is the MS running mass of the decoupling quark. 




o 1 1 1 1 1 1 

1 2 3 4 5 

Q , k x (GeV) 

Figure 2: QCD couplings in NNLO for a fixed number of active flavors rif = 3 and h.Q CD = 0.7 GeV. Dashed 
line (green): standard coupling a(Q 2 ); dotted line (blue): ghost-less or analytic coupling a(Q 2 ); continuous line (red): 
cascade or time-like coupling a(k 2 L ). 

Let us end this section summarizing the basic steps taken in the construction of the effective coupling for 
the gluon cascade of our model: 

1. subtraction of the Landau pole from the ordinary QCD coupling; 

2. inclusion of the absorptive effects related to the decay of time-like gluons in the coupling controlling jet 
evolution. 



3.1 Coupling in the DMW model 

Let us now evaluate the quantity 

1 f 1 ' 

a = — / dk ± &(k 2 ,) (48) 
m Jo 

parameterizing the leading non-perturbative effects in the Dokshitzer-Marchesini- Webber (DMW) model [S]- 
With fij = 2 GeV we obtain in our model «o = 0.40 for a(mz) = 0.12 and ao — 0.44 for a(mz) — 0.125. 
In general, we find that «o is roughly linear in a(mz)- A fit to e + e~ shape variables data using next-to-leading 
resummed formulas gives ao ~ 0.45 9 . 

9 G. Salam: private communication. 
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4 Threshold Resummation with Effective Coupling 



In this section we describe a model for threshold resummation in semi-inclusive beauty decays based on the 
effective coupling considered in the previous section. Basically, we replace in the resummation exponent the 
standard coupling with the effective one. 

4.1 iV-space 

In order to factorize multiple soft-gluon kinematic constraints, a transformation to ./V-space is required: 

a N (a) = [ (1 -t) N ~ l (T(t] a)dt, (49) 



JO 

where 

„(t; a) = S(t) + 0(a% (50) 

is the differential QCD form factor in the notation of 0. C F = {N^ - l)/(2N c ) = 4/3 with N c = 3 the 
number of colors and the plus distributions are defined as usual as: 

l 



P(t)+ = Pit) - S(t) / P(t')dt'. (51) 
Jo 

The form factor has an exponential form in iV-space: 

a N (a) = e G ^ a \ (52) 

where the exponent of the form factor reads: 

G N (a) = t ^ [(1 - yf- 1 1] { f " + S [5(Q 2 y)] + Z> [fi(QV)] ) , (53) 

with Q = wniB being the hard scale and w = 2Ex/ms- The functions A(<5), B{a) and D(a) have expansions 
in powers of the effective coupling: 



A(a) = ^2 A„ a n = A 1 a + A 2 a 2 + A 3 a 3 + ■ ■ ■ ; (54) 

n=l 

oo 

B(a) = ^B n Q" = Bia + B 2 Q 2 + '"; (55) 
n=l 

OO 

D(a) = A. «" = £>i a + D 2 a 2 + ■ ■ ■ . (56) 

n=l 

The resummation constants for the cascade coupling are obtained from the standard ones (usually in the MS 
scheme) by imposing equalities such as: 

A(a) = A(a), (57) 

where 

oo 

A(a) = Y A n a " = Aia + A 2 a 2 + A 3 a 3 H (58) 
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is the standard double- logarithmic function 10 . Expressing the cascade coupling in terms of the standard one, 
according to eq. (|44[l. we obtain: 11 

M = Ax; (59) 
A 2 = A 2 ; (60) 

A 3 = A 3 + ^l Al ; (61) 

M = A i + l (n/3 ) 2 A 2 + ~ ^ (tt/3o) 2 (62) 
3 6 fa 

The first two coefficients A\ and A 2 are the same for both couplings a and a, while the third-order one A3 
is modified going to the time-like coupling by a contribution proportional to the first-order coefficient. For 
rif = 3, A 3 k 1 is larger than A 3 « 0.3 by a factor 3, but it is still acceptably small. 

We have found that the inclusion of the NNLO terms involving A3, B 2 and D 2 — in particular A3 - is 
crucial for a good description of the experimental data. The NLO spectra are indeed peaked at too low hadron 
invariant masses and a sizable and positive value for A3 suppresses the elastic region and shifts the spectra to 
higher nix's, The model could be improved by including NNNLO terms, which require the knowledge of the 
coefficients A4, B3 and D3; at present, only B3 is analytically known |12| . 

Our model has been constructed by means of a power expansion in a single (effective) coupling a, i.e. higher 
orders are proportional to a™. In and in |13j a non-power expansion had been proposed involving a different 
coupling for any n, which has interesting theoretical properties. In second order (n = 2), for example, one has 
the coupling 

Pl{^ 2 + log 2 *i/A 2 ) 1 ; 

in place of a(/c^) 2 , with a(k\) given by eq. (|40|l . We have found that the non-power expansion does not offer 
a good description of the measured spectra. That is because 



a 2 {ki) -> while a(fci) 2 ~^ 7p ~ for k l ~> °- ( 64 ) 



That implies that second-order effects are suppressed in the soft region with the non-power expansion compared 
to the power expansion case. In general, the non-power expansion renders the higher-order effects very small 
P. But, as discussed above, in beauty decays, sizable third-order effects are needed to take the theoretical 
curves close to the data, disfavoring the non-power expansion. 

In order to include as many corrections as possible — higher order logiV terms, 1/N contributions, etc. - 
in agreement with the philosophy described in the introduction, we make the integration over y in Gn exactly, 
in numerical way. This is possible because the time- like coupling a(kj_), unlike the standard one, does not have 
the Landau singularity and is regular for any k\ > 0. 

4.2 Inverse Transform 

The form factor in momentum space is obtained by inverse transform: 

rC+ioo J AT 
C — too 



a) = / — {l-t)- N a N {a), (65) 



where the constant C is chosen so that the integration contour in the iV-plane lies to the right of all the 
singularities of <7Ar(a). In order to correctly implement multi-parton kinematics, the inverse transform from 

10 A compilation of the resummation constants in our normalization, with references to the original papers, can be found in 151 . 
^Analogous relations hold for the B^s and the Di's. 
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TV-space back to s-space is also made exactly in numerical way. Let us note that no prescription — such as 
the minimal prescription in the standard formalism 0] — is needed in our model because <7jy(a) is analytic for 
ReiV > 0. 



5 Radiative Decay 

The event fraction or partially-integrated rate for the radiative B decay 

B -> X s + 7 (66) 



can be written as [THUS] : 



where 



is a dimensionless variable, 



1 f* dT 



n u; Jo 



W L ~dY = K r( a )^& a ) + «) ( 67 ) 

(68) 



'_x_ 

> 2 B 



E(t; a) = I a(t'; a) dt' = 1 - %^ log 2 t - i ogt + (a 2 ) (69) 

ZTT 4 7T 



is the partially-integrated form factor and 

r (o) _ "mi Gp ml ml 2 
r ~ it 32tt3 ° 7 

is the lowest-order inclusive width, mi, ss 5 GeV is the beauty pole mass while m& is the MS mass evaluated in 
/i = mb- Their relation reads: 



m b = 



1 - aimb) ° F + O(qi) 



nib — 0.9 mft. (71) 



K r (a) is a short-distance coefficient function specific for this process and having an expansion in powers of a: 

K r {a) = 1 + aKl l) + a 2 K {2) + 0(a 3 ). (72) 
The explicit expression of the first-order term reads: 

where the Cj's are short-distance coefficient functions entering the effective b —> sj Hamiltonian, H. b ^ s ~ t , whose 
numerical values are given in the appendix, and the ?Vs are complex constants. D r (t; a) is a process-dependent 
remainder function, which is included to correctly describe also the high jet mass region t 0(1). In our 
leading-twist analysis, this function can be computed in perturbation theory and starts in first order: 

D r (t- a) = aD^(t) + a 2 £><?>(«) + 0(a 3 ), (74) 

with 

Di 1) (t) = lY:^fS\t). (75) 

i<j 7 

The /y(t)'s are functions whose explicit expression are given in the appendix. 
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Since the total width T r is infrared divergent beyond tree level, because of soft photon effects occurring in 
the spectrum for t —> 1, it is convenient to define an event fraction normalized to the partial rate 



(76) 



where S < 1 is a parameter. That is also convenient for experimental reasons: due to large backgrounds, the 
presently accessible range of hadron masses is at the most < t < 0.3 (see later). The event fraction normalized 
to T r (6) reads: 

*w-m£w«- (77) 

For S — ► 1, Rg(t) tends to the standard event fraction R r (t). The normalization condition is: 

Rs(S) = 1. (78) 
The differential spectrum is obtained by differentiation: 

A resummed expression of the following form holds: 

R s (t; a) = C s (a) E 4 (t; a) + D s (t; a), (80) 
where we have defined the form factor 

£(<5; a)' 



£,(f; a) = |^4, (81) 



which is normalized as 

X 5 (5; a) = 1. (82) 

The normalization condition (j78(l gives: 

C s (a) + D s {6; a) = 1. (83) 
The expansions of the coefficient function and the remainder function read: 

C s (a) = 1 + aCf ] + a 2 cf ] + 0(a 3 ); (84) 
D s (t;a) = aD ( l\t) + a 2 Df\t) + 0(a 3 ), (85) 

with 

^)-^- 

Note that we only expand D r (t\ a) in powers of a and not S(<5; a), because, for sufficient small 5, one can have 
a log 5 w O(l), implying need for resummation to any order in a. For 5 = 0.26, one obtains Cf ) ~ - 0.48, 
i.e. a 0(10%) correction to the treedevel coefficient function (E(t = 0.26; a) = 1.10). 



5.1 Phenomenology 

In Fig. |3| we compare the invariant hadron mass distribution for the radiative decay, dT r /dmx, obtained with 
our model with experimental data from the BaBar collaboration |16| . The data show a rather pronounced K* 
peak, which clearly cannot be accounted for in a perturbative QCD framework 12 . We have therefore discarded 

12 To have a point-to-point description of the data, one has to include by hand the contribution of this resonance, by means of 
one or more free parameters. 



13 



in the analysis the data points with mx < 1.1 GeV 13 . We obtain a minimum \ 2 = 12 for as(mz) = 0.1255 
for 13 data points, i.e. for 12 degrees of freedom (d.o.f.) because of the fixed normalization. Since we deal 
with the standard QCD coupling, let us write as from this section till the end of the paper. To improve the 
agreement of the theory with the data and to estimate the error on asijnz), we have performed a Gaussian 
smearing of Amx = 300 MeV of the data points and of the theoretical distribution, and we have discarded the 
points with mx < 800 MeV. We obtain a minimum % 2 = 6.8 for asijnz) = 0.1205 for 15 d.o.f.. By taking 
as an estimate of asijnz) the average of the above values and as an estimate of the error their difference, we 
quote: 14 

a s {m z ) = 0.123 ± 0.003 (m Xs ■ BABAR) . (87) 




In Fig. 01 we compare the photon energy spectrum computed in the framework of our model with data from 
the Cleo Collaboration ^7]. In the B rest-frame, 

t = 1 - (88) 
m B 

The photon energies are however measured in the T(4S') rest frame, in which the B mesons have a small, non- 
relativistic motion. In order to model the Doppler effect, we have convoluted the theoretical curve for - 
computed with a B meson at rest — with a normal distribution of a 1 = 150 MeV, as suggested by Cleo itself. 
Let us note that the Doppler effect is sufficient to completely wash out the K* peak. We obtain a minimum 
X 2 = 3.8 for as(mz) = 0.117 for 7 d.o.f.. Assuming complete independence of the experimental points, we 
allow the \ 2 to raise by one unit to estimate the error and we obtain: 

a s {m z ) = 0.117 ± 0.004 {E 1 : CLEO, ct 7 = 150 MeV). (89) 

To check the modeling above of the Doppler effect, we have used the following method. We have converted the 
mxs distribution by BaBar above to a photon spectrum in the B rest-frame and we have convoluted it with 



13 m K * = 892 MeV and Tk* = 51 MeV. 

14 We have taken ni(, = 4.8 GeV, rh c = 1.4 GeV and fh s = 0.3 GeV in the decoupling relations. In general, changing the 
MS masses in a reasonable range modifies the theoretical predictions in a negligible way. Increasing the MS masses is roughly 
equivalent to a slight increase of ajfinz). That is because, lowering the renormalization scale, the QCD coupling rises faster for a 
smaller number of active flavors. 
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a normal distribution with a variable cr 7 , obtaining the points (xi, y^(cr 7 ), a^). We have then minimized the 
quantity 

with respect to er 7 , where (xi, yi, <Ji) are the Cleo data 15 . We have found a minimum of H(a~ l ) for a 1 ~ 100 MeV, 
which gives similar results to the analysis with a 1 = 150 MeV: 

a s (m z ) = 0.118 ± 0.003 (£ 7 : CLEO, cr 7 = 100 MeV) , (91) 

with x 2 = 3.4. As intuitively expected, reducing ct 7 produces a shaper spectrum. 



CLEO » i 

8 - 



6 - 




-0.05 0.05 0.1 0.15 0.2 0.25 

t = m x 2 / m B 2 

Figure 4: B — > V s 7 photon spectrum from Cleo compared to our model. Dotted line (blue): as(mz) = 0.118 and 
o~ t — 100 MeV to model the Doppler effect (see text); continuous line (black): as(mz) = 0.117 and <j 7 = 150 MeV. 

In Fig. 03 we compare the predictions of our model with a spectrum from the BaBar collaboration jTHj. The 
BaBar spectrum is somewhat softer than the Cleo one — even though the difference is within one standard 
deviation; we have interpreted this difference as a resolution effect and we have convoluted our theoretical curve 
with a normal distribution with a slightly larger standard deviation, er 7 = 200 MeV 16 . We obtain a minimum 
X 2 = 5.1 for 9 d.o.f.. Performing a similar analysis as for the Cleo data, we obtain: 

a s {m z ) = 0.129 ± 0.005 {E 1 : BABAR, ct 7 = 200 MeV). (92) 

Following the minimization procedure above for the Cleo spectrum (see eq. (I90f) 'l. we obtain also for the BaBar 
photon spectrum a 1 ~ 100 MeV, to give: 

a s (m z ) = 0.130 ± 0.008 (£ 7 : BABAR, ct 7 = 100 MeV), (93) 

with x 2 — 8.0. Let us note that <r 7 and us(mz) are slightly anti-correlated because by increasing as{mz) 
more radiation is emitted with a smearing effect similar to the one of increasing ct 7 . 

The same analysis on the BaBar photon spectrum can be repeated for the Belle one (see Fig. EJ. The 
minimization of H{a 1 ) gives in this case <r 7 ~ 200 MeV. We obtain a minimum of x 2 — 5.3 for 8 d.o.f., to give 

a s (m z ) = 0.130 ± 0.005 {E 1 : BELLE, ct 7 = 200 MeV) . (94) 



Let us remark however that the two spectra entering eq. I9UI are independent on each other. 

A more sophisticated analysis from the experimentalists, including the true resolution functions, is strongly encouraged! 
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The over-all picture is that there is a good agreement of our model with the data in the region rax > 1 GeV, 
below which single resonances such as K and K* are expected to have a substantial effect in the dynamics. The 
extracted values of ots(mz) are in agreement with the world average 



a s (m z ) = 0.1176 ± 0.0020 
at most within two standard deviations. 



(PDG06) . 



(95) 



6 Semileptonic decay 

Resummed expressions for the triple-differential distribution in the inclusive charmless semileptonic B decays, 

B -> X u + I + v, (96) 

as well as for many double and single distributions have been given in ^] 03 E3 12] , so we do not repeat them 
here 17 . To compare with semileptonic data, we just supplement these resummed expressions with the QCD 
form factor a computed within our model. 

6.1 Hadron mass distribution 




Figure 7: 
0.119. 



1 .5 2 2.5 

m x (GeV) 

invariant hadron mass distribution in semileptonic decays from BaBar compared to our model for as(mz) 



In Fig. we compare the invariant hadron mass distribution in the semileptonic decay (|96[l in our model 
with data from the BaBar collaboration 23 . We discard the point with mx < 400 MeV, which is dominated 
by the ir peak, and the points with mx > 2.6 GeV, which give basically no information on the signal 18 . We 
obtain a minimum \ 2 — 1-1 f° r 5 d.o.f. and, using the method of the previous section, we obtain: 

as(mz) = 0.119 ± 0.003 (m Xu BABAR) . (97) 

Since the p width is larger than that of the K*, T p ~ 150 MeV ~ 3Tk*, and the binning is rather large 
(Amx = 400 MeV), we do not apply any smearing procedure in this case. 

17 A slightly different formalism, which uses light-cone variables and is equivalent to ours in leading twist, has been developed in 

m 

18 Let us stress that for mx > 1.7 GeV experimental errors become very large because of the large background coming from 
semileptonic b — » c transitions which have mx > mjj = 1.867 GeV. 
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Figure 
0.123. 



m x (GeV) 

invariant hadron mass distribution in semileptonic decay from Belle compared to our model for as(rnz) = 



In Fig. |H|we present a similar plot with Belle data |23|. To extract as(mz), we discard the first 7 points, 
having mx < 0.8 GeV. We obtain a minimum \ 2 — 5.3 with 7 d.o.f., to give as{mz) = 0.123 ± 0.006. Since 
the binning is smaller for Belle (Amj = 120 MeV) than for BaBar, the p peak is pretty visible now. To reduce 
the resonance effect, we convolve our theoretical curve and the experimental data with a normal distribution of 
a = 300 MeV, as we have made with the mx s spectrum in the previous section. By discarding the first four 
points, we obtain a minimum \ 2 — 0.41 for 10 d.o.f. to give as(rnz) = 0.115 ± 0.004. Combining the above 
measures as we have made for the mx s spectrum, we quote: 



a s {m z ) = 0.119 ± 0.004 



(mxu ■ BELLE) . 



Let us note that semileptonic distributions peak at smaller hadron masses than radiative ones because they 
have a smaller average hadron energy 9 : 



{E X ) si 0.7{E X )rd- 

For as{mz) — 0.123, we find for the peak positions in our model: 

mx„ f» 1.3 GeV while mx s ~ 1-7 GeV. 



(99) 



(100) 



We end this section by noting that the extracted values of as(mz) from the above measurements are in 
agreement with each other as well with the reference value within one standard deviation. 



6.2 Electron spectrum 



The electron spectrum in the decay l|9"tj|) is affected by a large background for 

™2 



E K < 



m B 



1 - 



2.31 GeV 



coming from the decays 



B -> X c + I + v. 



(101) 
(102) 
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This background is larger than the signal by two orders of magnitude because |K<fe| 2 /|K;6| 2 ~ 1CP 2 . In order to 
avoid the large errors coming from its subtraction, we have normalized the theory and the data to one in the 
region E e > 2.31 GeV. Instead of the electron energy, we prefer to use the variable 

x e = l -, 103 

m B 

which is equal to zero for the largest electron energy. The charm background occurs for x e > 0.125. To include 
the Doppler effect, we convolve our spectra with a normal distribution of standard deviation a e = 100 MeV. 




0.1 0.15 

Xg = 1 - 2 E e / m B 

Figure 9: electron spectrum in semileptonic charmless B decay from Cleo compared to our model with as(mz) = 0.118. 
The data and the theory are normalized to one in the charm background free region < x e < 0.125. 

In Fig. ED we compare the electron spectrum in our model 19 with data from the Cleo Collaboration j2S]- We 
obtain a minimum \ 2 — 30 for 13 d.o.f. and with the analysis described above we obtain: 



a s (m z } 



0.117 ± 0.005 



(E e : CLEO) 



(104) 



The over-all agreement of the model with the data is acceptable in all the measured range of electron energies. 
In the region affected by the charm background, experimental errors become however very large. 

In Fig. 1101 we compare our prediction with the electron spectrum measured by the BaBar collaboration 
|26| . In the \ 2 analysis we remove the 4 points with the smallest electron energies, which are affected by the 
subtraction of the charm background. We obtain a minimum \ 2 — 16 for 9 d.o.f. and we obtain 



a s {m z ) 



0.119 ± 0.005 



[E e : BABAR) 



(105) 



In Fig. ^] we compare our model with Belle data For the % 2 analysis we have discarded the seven 

points with the largest x e , i.e. with smallest electron energies. We obtain a minimum x 2 — 7 for 8 d.o.f. for 
ots{iTiz) ~ 0.135. Since the \ 2 is a rather irregular function of as(mz) in this case 
minimum — we are not able to estimate the error. 



without a well-shaped 



Finally, in Fig. 1121 we compare our model for as(mz) = 0.119 with a preliminary measure of the electron 
spectrum of the BaBar collaboration extending down to E e = 1.1 GeV [2E1- As it is clearly seen, the theoretical 
spectrum is harder than the experimental one. We do not known whether this discrepancy is related to a 
deficiency of our model or to an under-estimated charm background. 



19 Let us note that the tree-level electron spectrum has a maximum for x e = 0, at the largest electron energy, where it is flat. 
The shift of the maximum inside the kinematical domain, in x e ss 0.2, is a Sudakov effect. Because of infrared divergencies, soft 
radiation is always emitted and the high energy electron recoils against a neutrino and a massive up-quark jet, instead of a massless 
one. 
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0.1 0.15 

Xg = 1 - 2 E e / m B 

Figure 10: electron spectrum in semileptonic decay from BaBar compared to our model with as(mz) 
data and the theory are normalized to one in the charm background free region < x e < 0.125. 



0.119. The 




Figure 11: electron spectrum in semileptonic decay from Belle compared to our model for as(mz) 
and the theory are normalized to one in the charm background free region < x e < 0.125. 



0.135. The data 
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BABAR » i 




Xg = 1 - 2 E e / m B 

Figure 12: electron spectrum in semileptonic decay extending down to 1.1 GeV from BaBar (preliminary) compared to 
our model with as(mz) = 0.119. The data and the theory are normalized to one in the charm background free region 
< x e < 0.125. 

We may summarize our analysis of the electron spectra by saying that the agreement theory-data is less clear 
in this case. The agreement is acceptable in the charm background free region, i.e. for 2.31 < E e < 2.64 GeV; 
the errors in the measure of as(mz) are however larger than in previous cases. There is not a good agreement 
with the preliminary BaBar spectrum for small electron energies: our model predicts a broad maximum around 
E e — 2.1 GeV, while the data seem to peak at lower energies. 



7 Conclusions 

We have presented a model for the QCD form factor describing radiative and semileptonic B decay spectra 
based on soft-gluon resummation to next-to-next-to-leading logarithmic accuracy and on a power expansion 
in an ghost-less time-like coupling. The latter is free from infrared singularities (Landau ghost) and resumes 
absorptive effects in gluon cascades to all orders. 

The agreement with invariant hadron mass distributions in radiative and semi-leptonic decays measured by 
Cleo, BaBar and Belle is in general a good one. The x 2 /d.o.f. values are acceptable and the extracted values 
of ats(mz) are in agreement with the current PDG average within two at most standard deviations. 

The agreement with the electron spectra in semi-leptonic decays is in general is, in general, not as good. 
Even restricting the analysis to the end-point region free from the charm background (2.31-^2.64 GeV), x 2 /d.o.f. 
values are larger and the extracted values of as(mz) are generally less accurate than in previous cases. The 
preliminary BaBar measure of the electron spectrum down to 1.1 GeV is not in good agreement with our 
model, which predicts a harder spectrum, with a broad maximum around 2.1 GeV. We do not know whether 
the discrepancy is to be attributed to a deficiency of our model or to an under-subtracted charm background. 
In the former case, one could think to a non-perturbative component which is accidentally larger in the electron 
spectrum than in other semi-leptonic or radiative spectra. 

In general, the model seems to work quite well, validating the idea that Fermi-motion effects can be described 
in a resummed pQCD framework with an effective QCD coupling, which remains reasonably smaller than one 
in the relevant integration domain. Since the effective coupling is constructed by means of an extrapolation at 
low energy of the standard coupling, non-perturbative Fermi-motion effects are connected in a smooth way to 
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the perturbative ones — namely soft gluon radiation — in our model. More accurate data on any distribution 
sensitive to soft- gluon effects are needed to put the model to a stringent test. Theoretical predictions could 
be sharpened in the future by including second-order corrections to the coefficient functions and remainder 
functions, as soon as they become available; that would allow to work within a complete NNLO approximation. 

We have found that the inclusion of the NNLO effects in our model is crucial for a good description of 
the data; the model could be improved by including NNNLO terms. We have also found that the non-power 
expansion proposed in P3E] does not accurately describe soft-gluon effects. 

Let us end with a general comment. We find it remarkable that with such a simple model as the one we 
have formulated, without any adjustable parameter, it is possible to extract reasonable values of as{mz) from 
spectra with a hard scale of just a few GeV. 
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A Coefficient Functions and Remainder Functions for the Radiative 
Decay 

The coefficient functions C, of the local operators Oi entering the effective Hamiltonian 

8 

n b-> sl =J2Ci Oi (106) 

i=i 

can be taken as: 

(Ci, C 2 , C 3 , C 4 , C 5 , C 6 , C 7 , C 8 ) = (-0.480, 1.023, -0.0045, -0.0640, 0.0004, 0.0009, -0.304, -0.148). 

(107) 

The analytic expressions of the C^'s as functions of nib, m t, and asijnz) as well as of the O^s can be found 
in [221 20 ■ The coefficients Ci for i — 3,4,5,6 are very small, implying that the contributions of the related 
operators can be neglected. 

The functions entering the leading-order remainder function D^\t) read in our conventions: 



hi(t) = +_/ 22 (t); (108) 

36 

/«(*) = -i/aa(*)i (109) 

fir® = - 1 /ar(*)5 (HO) 

fm(t) = --/«.(*); (HI) 
6 

Unlike previous works ( 1141 and 1151 '). we always insert in the formulas the corrected C7 = + os/(47r) . 
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Mt) = 

Mt) = 

Mt) = 

frr(t) = 

/78(f) = 

/«»(*) = 
We have defined: 
and 
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27 



k i t 



(i-t)/fc 



dv (1 — kv) 







G(v) 1 
w 2 



2 M/fc 



+ / - fc«r 



9 1 



-g/a7(*)5 



(l-i)A 



(iu Re 



G(v) 



1/k 



(l-t)/fc 



(l-t)/fc 



efe(l - few) Re |G(<u) + - 



G(«) 


1 


1 


V 


+ 2 





- [30 + 3i-2t 2 -3 (4-t) logt] ; 
9 

i [2tt 2 -27i + 3i 2 -i 3 + 12tlogt-12Li 2 (l-f)] 



(112) 

(113) 

(114) 
(115) 
(116) 



if - 27r 2 + i [21 + (9 -2*) t] + 24 log(l - 1) - 6 log — [t (2 + t) + A logfl - f)] + 
81 L m s 

- 3* (2 + *) logi + 12Li 2 (l-i)|. (117) 



G( V ) 



2 arctan 2 

2 



2 log 



2 



27TZ log 



\/v + \/v — 4 
2 



v < 4: 



n, — - ri . 

mf 



(118) 



(119) 



B QCD form factor 



In this appendix we tabulate the values of the QCD form factor a(u; w) in our model as a function of the 
infrared variable 



■ - x 



E X + JW X - m x *E% 
and of the total final hadron energy 



(for m x < E x ) (120) 



w ee ^ (121) 

TUB 

for as(mz) = 0.115, 0.120 and 0.125. The hard scale in the process is 

Q = wm B . (122) 

In the radiative case one sets t = u and w = 1 while in the semileptonic case the form factor as a function of w 
(0 < w < 1) is needed |9l 1201 |2T] , By using the following tables, the reader can obtain the form factor for all 
the values of u and w by means of a straightforward interpolation, avoiding the delicate numerical integrations 
related to the Mellin transform and to the inverse Mellin transform. In agreement with physical intuition, by 
lowering the hard scale, the peak of the form factor broadens and shifts to larger u's because of the coupling 
growth. 
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<r(u; w) : a s (m z ) = 0.115 



u 


w = 0.10 


w = 0.28 


w = 0.46 


w = 0.64 


w = 0.82 


w = 1.00 


0.01 


2.89 x 10" b 


5.094 x 10" b 


1.218 x 10" b 


1.088 x 10" 4 


8.073 x 10" 4 


3.825 x 10" 3 


0.02 


3.063 x 10- e 


1.028 x 10" 4 


3.525 x 10" 3 


3.077 x 10" 2 


1.344 x 10" 1 


3.923 x 10" 1 


0.03 


3.686 x 10" b 


3.574 x 10" 3 


6.732 x 10" 2 


3.688 x 10" 1 


1.12 


2.426 


0.04 


4.715 x 10" 5 


3.068 x 10" 2 


3.597 x 10" 1 


1.407 


3.295 


5.788 


0.05 


4.848 x 10" 4 


1.277 x 10" 1 


1.03 


3.102 


5.948 


8.903 


0.06 


2.577 x 10" 3 


3.485 x 10" 1 


2.065 


5.029 


8.209 


1.082 x 10 1 


0.07 


9.099 x 10" 3 


7.256 x 10" 1 


3.313 


6.749 


9.64 


1.145 x 10 1 


0.08 


2.451 x 10" 2 


1.256 


4.578 


8.002 


1.021 x 10 1 


1.111 x 10 1 


0.09 


5.449 x 10" 2 


1.905 


5.697 


8.713 


1.008 x 10 1 


1.019 x 10 1 


0.10 


1.052 x 10" 1 


2.621 


6.568 


8.929 


9.493 


9.004 


0.11 


1.821 x 10" 1 


3.346 


7.152 


8.754 


8.637 


7.748 


0.12 


2.894 x 10" 1 


4.031 


7.456 


8.304 


7.661 


6.549 


0.13 


4.296 x 10" 1 


4.636 


7.511 


7.681 


6.669 


5.469 


0.14 


6.029 x 10" 1 


5.135 


7.366 


6.964 


5.725 


4.531 


0.15 


8.075 x 10" 1 


5.516 


7.067 


6.213 


4.863 


3.734 


0.16 


1.04 


5.776 


6.661 


5.472 


4.099 


3.065 


0.17 


1.295 


5.92 


6.186 


4.769 


3.435 


2.51 


0.18 


1.568 


5.958 


5.673 


4.119 


2.864 


2.05 


0.19 


1.851 


5.902 


5.145 


3.532 


2.379 


1.67 


0.20 


2.138 


5.769 


4.621 


3.008 


1.968 


1.356 


0.21 


2.423 


5.571 


4.114 


2.548 


1.621 


1.097 


0.22 


2.7 


5.325 


3.635 


2.146 


1.329 


8.814 x 10" 1 


0.23 


2.963 


5.042 


3.187 


1.798 


1.083 


7.027 x 10" 1 


0.24 


3.208 


4.734 


2.776 


1.498 


8.765 x 10" 1 


5.537 x 10" 1 


0.25 


3.43 


4.411 


2.401 


1.24 


7.021 x 10" 1 


4.293 x 10" 1 


0.26 


3.627 


4.081 


2.063 


1.019 


5.55 x 10" 1 


3.249 x 10" 1 


0.27 


3.797 


3.75 


1.761 


8.29 x 10" 1 


4.307 x 10" 1 


2.372 x 10" 1 


0.28 


3.937 


3.423 


1.491 


6.665 x 10" 1 


3.255 x 10" 1 


1.632 x 10" 1 


0.29 


4.048 


3.105 


1.253 


5.273 x 10" 1 


2.363 x 10" 1 


1.005 x 10" 1 


0.30 


4.128 


2.798 


1.042 


4.079 x 10" 1 


1.606 x 10" 1 


4.656 x 10" 2 


0.31 


4.18 


2.504 


8.571 x 10" 1 


3.057 x 10" 1 


9.62 x 10" 2 


1.027 x 10" 4 


0.32 


4.202 


2.227 


6.948 x 10" 1 


2.18 x 10" 1 


4.139 x 10" 2 


-3.977 x 10" 2 


0.33 


4.197 


1.965 


5.528 x 10" 1 


1.427 x 10" 1 


-5.313 x 10" 3 


-7.377 x 10" 2 


0.34 


4.166 


1.722 


4.29 x 10" 1 


7.814 x 10" 2 


-4.524 x 10" 2 


-1.026 x 10" 1 


0.35 


4.112 


1.496 


3.211 x 10" 1 


2.272 x 10" 2 


-7.955 x 10" 2 


-1.27 x 10" 1 


0.36 


4.035 


1.287 


2.274 x 10" 1 


-2.485 x 10" 2 


-1.093 x 10" 1 


-1.476 x 10" 1 


0.37 


3.938 


1.096 


1.463 x 10" 1 


-6.568 x 10" 2 


-1.353 x 10" 1 


-1.649 x 10" 1 


0.38 


3.824 


9.216 x 10" 1 


7.613 x 10" 2 


-1.007 x 10" 1 


-1.578 x 10" 1 


-1.796 x 10" 1 


0.39 


3.693 


7.63 x 10" 1 


1.552 x 10" 2 


-1.307 x 10" 1 


-1.771 x 10" 1 


-1.918 x 10" 1 


0.40 


3.55 


6.195 x 10" 1 


-3.678 x 10" 2 


-1.564 x 10" 1 


-1.936 x 10" 1 


-2.021 x 10" 1 


0.50 


1.791 


-1.934 x 10" 1 


-2.85 x 10" 1 


-2.775 x 10" 1 


-2.596 x 10" 1 


-2.42 x 10" 1 


0.60 


3.035 x 10" 1 


-3.989 x 10" 1 


-3.229 x 10" 1 


-2.851 x 10" 1 


-2.551 x 10" 1 


-2.346 x 10" 1 


0.70 


-4.713 x 10" 1 


-4.072 x 10" 1 


-3.088 x 10" 1 


-2.615 x 10" 1 


-2.335 x 10" 1 


-2.15 x 10" 1 


0.80 


-7.243 x 10" 1 


-3.596 x 10" 1 


-2.718 x 10" 1 


-2.311 x 10" 1 


-2.076 x 10" 1 


-1.92 x 10" 1 


0.90 


-6.43 x 10" 1 


-2.953 x 10" 1 


-2.286 x 10" 1 


-1.975 x 10" 1 


-1.792 x 10" 1 


-1.653 x 10" 1 


0.98 


-4.249 x 10" 1 


-2.458 x 10" 1 


-2.017 x 10" 1 


-1.803 x 10" 1 


-1.672 x 10" 1 


-1.58 x 10" 1 
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<r(u; w) : a s (m z ) = 0.120 



u 


w = 0.10 


w = 0.28 


w = 0.46 


w = 0.64 


w = 0.82 


w = 1.00 


0.01 


2.575 x 10" b 


4.491 x 10" b 


6.691 x 10" b 


2.153 x 10" b 


1.321 x 10" 4 


6.542 x 10" 4 


0.02 


2.748 x 10" b 


2.286 x 10" 5 


7.371 x 10" 4 


7.231 x 10" 3 


3.539 x 10" 2 


1.147 x 10" 1 


0.03 


2.992 x 10" e 


8.639 x 10" 4 


1.881 x 10" 2 


1.185 x 10" 1 


4.081 x 10" 1 


9.901 x 10" 1 


0.04 


1.341 x 10" 5 


8.992 x 10" 3 


1.247 x 10" 1 


5.674 x 10" 1 


1.517 


2.995 


0.05 


1.359 x 10" 4 


4.372 x 10" 2 


4.23 x 10" 1 


1.493 


3.282 


5.527 


0.06 


8.021 x 10" 4 


1.357 x 10" 1 


9.753 x 10" 1 


2.797 


5.245 


7.775 


0.07 


3.121 x 10" 3 


3.155 x 10" 1 


1.761 


4.239 


6.961 


9.285 


0.08 


9.146 x 10" 3 


6.013 x 10" 1 


2.695 


5.582 


8.182 


9.98 


0.09 


2.192 x 10" 2 


9.932 x 10" 1 


3.67 


6.663 


8.849 


9.991 


0.10 


4.521 x 10" 2 


1.475 


4.587 


7.407 


9.021 


9.516 


0.11 


8.315 x 10" 2 


2.018 


5.373 


7.812 


8.805 


8.74 


0.12 


1.397 x 10" 1 


2.589 


5.986 


7.916 


8.316 


7.817 


0.13 


2.181 x 10" 1 


3.157 


6.41 


7.774 


7.657 


6.855 


0.14 


3.208 x 10" 1 


3.692 


6.65 


7.447 


6.911 


5.924 


0.15 


4.489 x 10" 1 


4.171 


6.724 


6.991 


6.141 


5.062 


0.16 


6.024 x 10" 1 


4.58 


6.657 


6.453 


5.388 


4.29 


0.17 


7.799 x 10" 1 


4.909 


6.474 


5.873 


4.681 


3.612 


0.18 


9.791 x 10" 1 


5.155 


6.201 


5.282 


4.033 


3.024 


0.19 


1.197 


5.317 


5.862 


4.702 


3.451 


2.521 


0.20 


1.429 


5.401 


5.478 


4.149 


2.935 


2.093 


0.21 


1.671 


5.413 


5.067 


3.633 


2.482 


1.73 


0.22 


1.918 


5.361 


4.643 


3.159 


2.089 


1.423 


0.23 


2.167 


5.254 


4.219 


2.729 


1.748 


1.163 


0.24 


2.412 


5.1 


3.803 


2.343 


1.454 


9.444 x 10" 1 


0.25 


2.649 


4.909 


3.404 


1.999 


1.201 


7.596 x 10" 1 


0.26 


2.875 


4.689 


3.025 


1.695 


9.845 x 10" 1 


6.034 x 10" 1 


0.27 


3.085 


4.446 


2.67 


1.427 


7.986 x 10" 1 


4.712 x 10" 1 


0.28 


3.278 


4.186 


2.34 


1.193 


6.394 x 10" 1 


3.593 x 10" 1 


0.29 


3.451 


3.916 


2.037 


9.873 x 10" 1 


5.031 x 10" 1 


2.641 x 10" 1 


0.30 


3.602 


3.64 


1.761 


8.082 x 10" 1 


3.864 x 10" 1 


1.828 x 10" 1 


0.31 


3.73 


3.362 


1.51 


6.522 x 10" 1 


2.864 x 10" 1 


1.127 x 10" 1 


0.32 


3.834 


3.086 


1.283 


5.165 x 10" 1 


2.008 x 10" 1 


5.236 x 10" 2 


0.33 


3.913 


2.815 


1.08 


3.986 x 10" 1 


1.274 x 10" 1 


5.426 x 10" 4 


0.34 


3.968 


2.552 


8.978 x 10" 1 


2.963 x 10" 1 


6.449 x 10" 2 


-4.377 x 10" 2 


0.35 


3.999 


2.298 


7.355 x 10" 1 


2.076 x 10" 1 


1.037 x 10" 2 


-8.151 x 10" 2 


0.36 


4.007 


2.055 


5.915 x 10" 1 


1.308 x 10" 1 


-3.641 x 10" 2 


-1.136 x 10" 1 


0.37 


3.992 


1.824 


4.64 x 10" 1 


6.426 x 10" 2 


-7.71 x 10" 2 


-1.407 x 10" 1 


0.38 


3.956 


1.607 


3.515 x 10" 1 


6.723 x 10" 3 


-1.125 x 10" 1 


-1.636 x 10" 1 


0.39 


3.9 


1.403 


2.527 x 10" 1 


-4.299 x 10" 2 


-1.433 x 10" 1 


-1.83 x 10" 1 


0.40 


3.826 


1.212 


1.659 x 10" 1 


-8.59 x 10" 2 


-1.697 x 10" 1 


-1.993 x 10" 1 


0.50 


2.42 


-9.317 x 10" 3 


-2.727 x 10" 1 


-2.938 x 10" 1 


-2.836 x 10" 1 


-2.655 x 10" 1 


0.60 


7.765 x 10" 1 


-4.072 x 10" 1 


-3.595 x 10" 1 


-3.213 x 10" 1 


-2.858 x 10" 1 


-2.61 x 10" 1 


0.70 


-3.087 x 10" 1 


-4.693 x 10" 1 


-3.545 x 10" 1 


-2.965 x 10" 1 


-2.617 x 10" 1 


-2.388 x 10" 1 


0.80 


-7.857 x 10" 1 


-4.228 x 10" 1 


-3.12 x 10" 1 


-2.604 x 10" 1 


-2.312 x 10" 1 


-2.12 x 10" 1 


0.90 


-7.863 x 10" 1 


-3.432 x 10" 1 


-2.588 x 10" 1 


-2.199 x 10" 1 


-1.975 x 10" 1 


-1.808 x 10" 1 


0.98 


-5.003 x 10" 1 


-2.753 x 10" 1 


-2.213 x 10" 1 


-1.96 x 10" 1 


-1.805 x 10" 1 


-1.699 x 10" 1 
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<r(u; w) : a s (m z ) = 0.125 



u 


w = 0.10 


w = 0.28 


w = 0.46 


w = 0.64 


w = 0.82 


w = 1.00 


0.01 


2.309 x 10" b 


4.039 x 10" b 


5.47 x 10" b 


8.764 x 10" b 


2.845 x 10" b 


1.247 x 10" 4 


0.02 


2.46 x 10" b 


8.166 x 10" b 


1.706 x 10" 4 


1.785 x 10" 3 


9.505 x 10" 3 


3.339 x 10" 2 


0.03 


2.641 x 10" b 


2.309 x 10" 4 


5.49 x 10" 3 


3.839 x 10" 2 


1.458 x 10" 1 


3.871 x 10" 1 


0.04 


5.702 x 10" b 


2.809 x 10" 3 


4.39 x 10" 2 


2.243 x 10" 1 


6.666 x 10" 1 


1.448 


0.05 


4.474 x 10" b 


1.562 x 10" 2 


1.727 x 10" 1 


6.899 x 10" 1 


1.695 


3.153 


0.06 


2.796 x 10" 4 


5.422 x 10" 2 


4.501 x 10" 1 


1.47 


3.092 


5.072 


0.07 


1.174 x 10" 3 


1.387 x 10" 1 


9.02 x 10" 1 


2.485 


4.589 


6.778 


0.08 


3.689 x 10" 3 


2.875 x 10" 1 


1.512 


3.596 


5.938 


8.023 


0.09 


9.412 x 10" 3 


5.117 x 10" 1 


2.231 


4.665 


6.988 


8.735 


0.10 


2.055 x 10" 2 


8.124 x 10" 1 


2.996 


5.587 


7.678 


8.96 


0.11 


3.98 x 10" 2 


1.181 


3.745 


6.303 


8.017 


8.791 


0.12 


7.011 x 10" 2 


1.603 


4.429 


6.79 


8.05 


8.34 


0.13 


1.144 x 10" 1 


2.057 


5.012 


7.056 


7.838 


7.71 


0.14 


1.752 x 10" 1 


2.524 


5.473 


7.121 


7.447 


6.985 


0.15 


2.547 x 10" 1 


2.982 


5.805 


7.019 


6.935 


6.228 


0.16 


3.542 x 10" 1 


3.415 


6.01 


6.781 


6.354 


5.482 


0.17 


4.742 x 10" 1 


3.809 


6.098 


6.443 


5.743 


4.777 


0.18 


6.145 x 10" 1 


4.152 


6.082 


6.034 


5.133 


4.127 


0.19 


7.741 x 10" 1 


4.439 


5.975 


5.582 


4.544 


3.54 


0.20 


9.512 x 10" 1 


4.666 


5.794 


5.108 


3.99 


3.017 


0.21 


1.143 


4.832 


5.554 


4.63 


3.477 


2.558 


0.22 


1.348 


4.939 


5.268 


4.161 


3.011 


2.156 


0.23 


1.561 


4.99 


4.95 


3.71 


2.591 


1.808 


0.24 


1.781 


4.989 


4.611 


3.285 


2.216 


1.507 


0.25 


2.002 


4.942 


4.26 


2.889 


1.883 


1.247 


0.26 


2.223 


4.855 


3.907 


2.525 


1.59 


1.024 


0.27 


2.438 


4.731 


3.558 


2.192 


1.333 


8.328 x 10" 1 


0.28 


2.647 


4.577 


3.217 


1.891 


1.108 


6.686 x 10" 1 


0.29 


2.844 


4.399 


2.89 


1.619 


9.118 x 10" 1 


5.278 x 10" 1 


0.30 


3.029 


4.199 


2.579 


1.376 


7.41 x 10" 1 


4.067 x 10" 1 


0.31 


3.199 


3.984 


2.286 


1.159 


5.926 x 10" 1 


3.022 x 10" 1 


0.32 


3.352 


3.757 


2.012 


9.665 x 10" 1 


4.639 x 10" 1 


2.117 x 10" 1 


0.33 


3.486 


3.522 


1.758 


7.956 x 10" 1 


3.523 x 10" 1 


1.334 x 10" 1 


0.34 


3.601 


3.282 


1.523 


6.446 x 10" 1 


2.557 x 10" 1 


6.584 x 10" 2 


0.35 


3.696 


3.04 


1.308 


5.114 x 10" 1 


1.719 x 10" 1 


7.763 x 10" 3 


0.36 


3.769 


2.8 


1.112 


3.943 x 10" 1 


9.927 x 10" 2 


-4.2 x 10" 2 


0.37 


3.822 


2.563 


9.342 x 10" 1 


2.914 x 10" 1 


3.593 x 10" 2 


-8.451 x 10" 2 


0.38 


3.855 


2.331 


7.732 x 10" 1 


2.013 x 10" 1 


-1.942 x 10" 2 


-1.207 x 10" 1 


0.39 


3.866 


2.106 


6.284 x 10" 1 


1.224 x 10" 1 


-6.775 x 10" 2 


-1.515 x 10" 1 


0.40 


3.858 


1.89 


4.984 x 10" 1 


5.351 x 10" 2 


-1.098 x 10" 1 


-1.775 x 10" 1 


0.50 


2.917 


2.99 x 10" 1 


-2.15 x 10" 1 


-2.945 x 10" 1 


-3.031 x 10" 1 


-2.883 x 10" 1 


0.60 


1.277 


-3.637 x 10" 1 


-3.894 x 10" 1 


-3.592 x 10" 1 


-3.194 x 10" 1 


-2.902 x 10" 1 


0.70 


-6.465 x 10" 2 


-5.265 x 10" 1 


-4.05 x 10" 1 


-3.365 x 10" 1 


-2.936 x 10" 1 


-2.656 x 10" 1 


0.80 


-7.947 x 10" 1 


-4.951 x 10" 1 


-3.589 x 10" 1 


-2.94 x 10" 1 


-2.578 x 10" 1 


-2.344 x 10" 1 


0.90 


-9.333 x 10" 1 


-3.996 x 10" 1 


-2.936 x 10" 1 


-2.452 x 10" 1 


-2.179 x 10" 1 


-1.978 x 10" 1 


0.98 


-5.886 x 10" 1 


-3.088 x 10" 1 


-2.428 x 10" 1 


-2.128 x 10" 1 


-1.947 x 10" 1 


-1.826 x 10" 1 
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